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ABSTRACT 



We calculate gamma-ray spectra of pulsars for a set of electron acceleration models 
within the polar-cap scenario. We discuss our results in the light of the recent analysis 
(Kanbach 1999) of temporal properties of the Vela pulsar in the EGRET energy range. 

<N 

' 1. Introduction 

(N 

\ In a recent analysis of gamma-ray double-peak pulses of the Crab, the Vela and Geminga, 

as detected by EGRET , Kanbach (1999) pointed to an intriguing possibility that for the Vela 
pulsar the phase separation (A peak ) between the two peaks may actually be energy dependent. 
The energy-averaged value of the separation is very large in all three objects, between 0.4 and 
■ 0.5 (Fierro, Michelson & Nolan 1998), and the effect by itself is of the order of a few percent or 

less. In the case of the Vela pulsar (the object with the best photon statistics), the plot of A peak 
against energy (Fig. 2 middle panel of Kanbach 1999) shows that A peak is decreasing by about 5% 



o 



o 

over 20 energy intervals covering the range between ~ 50 MeV and ~ 9 GeV. The scatter of the 
points seems to us, however, also consistent with the separation staying at a constant level of 0.43, 
provided that we reject the lowest and the highest energy interval. 

X' 

$_i ■ The problem raised by Kanbach is important by itself from a theoretical point of view, 



regardless of whether his finding becomes well established or not. With future high-sensitivity 
missions like GLAST, any firm empirical relation between the peak-separation A peak and the 
photon energy e may serve as a tool to verify some models of pulsar activity. 

The presence of two peaks in gamma-ray pulses with large (0.4 - 0.5) phase separation may 
be understood within a scenario of a single canonical polar cap (e.g. Daugherty & Harding 1996, 
Miyazaki & Takahara 1997). One need to assume, however, a nearly aligned rotator, i.e. a rotator 
where three characteristic angles are of the same order: a - the angle between spin axis and the 
magnetic moment ft, 6-y - the opening angle between a direction of the gamma-ray emission and /J, 
and C - the angle between Vt and the line of sight. For a canonical polar cap and instant electron 
acceleration, # 7 roughly equals 0.02 /y/P radians only (where P denotes a spin period). To avoid 
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uncomfortably small characteristic angles, Daugherty & Harding (1996) postulated that primary 
electrons come from more extended polar caps, and with the acceleration occuring at a height h 
of several neutron-star radii R ns . The latter assumption may be justified by a GR effect found by 
Muslimov & Tsygan (1992). 

The aim of this paper is to present general properties of A peak (e) obtained numerically (and 
to some degree semi-analytically) for four simplified versions of a polar-cap activity model. In 
Section 2 we outline the model. Section 3 describes the results and offers some explanation of the 
presented effects. Conclusions follow in Section 4. 



2. The Model 

We use a polar cap model with beam particles (primary electrons) distributed evenly along 
a hollow cone formed by the magnetic field lines from the outer rim of a canonical polar cap, i.e. 
with an opening angle #i n it = pc , where 9 pc ~ (2ir R ns / cP) 1 / 2 radians at the stellar surface level 
(h = 0). In one case (model C - see below) we assume #; n it = 26? pc The essential ingredients of 
the high-energy processes had been introduced by Daugherty & Harding (1982), with high-energy 
radiation due to curvature and synchrotron processes (CR and SR, respectively) induced by 
primary electrons accelerated to ultrarelativistic energies. The pulsar parameters are those of the 
Vela: the spin period P = 0.0893 s, and the dipolar magnetic field at the polar cap B pc rs 10 12 G. 

Within the hollow-cone geometry we have considered three scenarios for electron acceleration: 
model A - beam particles are injected at a height h- m i t = with some initial ultrarelativistic energy 
-Emit (the values of £7 m it are listed in Table 1) and no subsequent acceleration; 
model B - similar to model A but b earn particles are injected at a height /^init — 1 Rns] 
model C - beam particles are injected at a height h- m \t = 2R ns with a low energy E{ n n and then 
accelerated by a longitudinal electric field £ present over a characteristic scale height Ah = 0.6i? ns , 
resulting in total potential drop Vq\ 

s(h) = i y °/ A/l ' for hinit - h - ^ init + A/l ) m 

[ ' \0, elsewhere. 1 ' 

For comparison, we considered a model with a uniform electron distribution over the entire 
polar cap surface (i.e. 9- m \t £ [0, pc ]): 

model D - beam particles are injected at a height h^ = with initial ultrarelativistic energy E{ n i t 
as in the model A, and no subsequent acceleration. 



The values of E[ n i t in models A, B and D, and the potential drop Vq in model C were chosen 
to yield similar number of secondary pairs - about 10 3 per beam particle. Table 1. summarizes 
the properties of the models. 
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Table 1. Model parameters. 





-Bpc 
[10 12 G] 


a 
[deg] 


/linit 
[Rns] 


dinit 


primary electrons 


Model A 


1.0 


3.0 


0.0 


1.0 


-Emit = 8.68 TeV, no acceleration 


Model B 


1.0 


5.0 


1.0 


1.0 


Einit = 20.0 TeV, no acceleration 


Model C 


3.0 


10.0 


2.0 


2.0 


-Emit = 0.5 MeV, acceleration (see eq.(l)) 
with Vo = 2.5 x 10 7 volts 


Model D 


1.0 


3.0 


0.0 


[0,1] 


Einit = 8.68 TeV, no acceleration 



3. Results 

We have calculated numerically the pulse shapes as a function of photon energy for all four 
models. The main difference between models A and B is due to different values of /ii n ;t which 
result in different locations of origin of secondary particles. Changing these locations is an easy 
way to modify spectral properties of emergent radiation and enables to change (preferrably - to 
increase) the angle a as constrained by the observed energy-averaged peak separation A pcak 0.43 
(Kanbach 1999). For a given model there are two possible values of the angle ( resulting in a 
desired peak separation A pea k at some energy e (which we choose to be ~ 1 GeV), and we took the 
larger one in each case: £ = 3.75,4.5, 15., and 3.65 degrees for models A, B, C, and D, respectively. 

In Fig. 1 we present dependency of peak separation on photon energy for all four models. 
Within a 'low-energy' range (below a few GeV) the peak separation either remains constant 
(Model B) or slightly decreases with the increasing photon energy (Models A, C, and D). At 
a critical energy e tU rn around a few GeV, the separation A peak undergoes a sudden turn: for 
e > etum it either rapidly increases (models A, B and C) or rapidly decreases (model D). For 
comparison, the overall trend (taking into account a substantial scatter of points) resulting from 
the analysis by Kanbach (see Fig.2, middle panel, of Kanbach 1999) is marked schematically. 

It would be easy to understand why A peak should stay at a constant level in the range e < e turn 
if the gamma-ray photons were due exclusively to CR (especially in models with no acceleration). 
A preferred direction of this emission would be set up already at the lowest magnetospheric 
altitudes. This is because a monotonic decrease of the electron's energy E, and initial increase in 
the dipolar curvature radius p cr make the level of contribution to the one-particle CR-spectrum 
at all energies to be the highest one just at the initial altitude h- m it- However, in the energy 
range below a few GeV the gamma-ray emission in our models is dominated by the synchrotron 
radiation (SR) due to secondary pairs in the cascades. The resulting behaviour of A pcak - 
either its slight decrease with increasing e (models A and C) or no change at all (model B) - is 
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Fig. 1. — Upper panels: Peak separation A p against photon energy e, as predicted by models A, 
B, C and D (red dots). Steep blue lines in the three left panels (hollow-cone models) are due to semi- 
analytical calculations for A peak (e) valid for e > eturn (see text). The overall observational trend, 
taking into account a substantial scatter of points (see Fig. 2, middle panel, of Kanbach 1999), is 
marked schematically for reference in all four panels as a green parallelogram. The assumed values 
of £ - the angle between fi and the line of sight - are also indicated. 

Lower panels: Energy output per logarithmic energy bandwidth at the first peak as a function of 
photon energy e for models A, B, C and D. [Note: these are not instantaneous spectra.] Units on 
the vertical axis are arbitrary. 
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actually mediated by several factors which influence directional and spectral properties of the SR. 
These include energy and pitch angle distributions of secondary pairs, as well as the spatial 
spread of these pairs within the magnetospere (see e.g. Rudak & Dyks 1999). These factors do 
change from one model to another. 

At the energy eturn the separation A peak undergoes a sudden turn and starts changing rapidly 
as e increases. (This occurs above ~ 1 GeV, where most photons are due to CR.) This is the 
regime where the position of each peak in the pulse is determined by the magnetospheric opacity 
due to jB — > e ± . 

For the hollow-cone models (A, B and C) the photons in both peaks of a pulse come from 
low magnetospheric altitudes with narrow opening angles. When e is high enough these photons 
will be absorbed by the magnetic field with subsequent pair-creation. In other words, inner parts 
of the 'original' peaks in the pulse will be eaten-up and the gap between the peaks (i.e. the 
peak separation) will increase. Photons which now found themselves in the 'new' peaks come 
from higher altitudes (the magnetosphere is transparent to them) and have wider opening angles. 
To quantify this line of arguments we present a simple semi-analytical solution for A peak as a 
function of e, which reproduces with astonishing accuracy our Monte Carlo results (see Fig.l). For 
each point (at a given radial coordinate r) lying on a magnetic field line with a known angular 
coordinate 9 on the stellar surface, one defines the "escape energy" (Harding et al. 1997), which 
is the upper limit e eS c for photon energy if the photon is to avoid magnetic absorption when 
propagating outwards. We have approximated e esc with a power-law formula £csc(0 — Q(^"/-Rns) 
MeV, and the values of a and b were found by fits to numerical solutions for each model. We 
found a ~ 7.83 x 10 2 , b ~ 2.49 for B = 10 12 G and 6 = 8 pc (models A and B); a ~ 1.25 x 10 2 , 
b ~ 2.50 for B = 3 x 10 12 G and = 20 pc (model C). Photons of some energy e will escape the 
magnetosphere only when they are emitted at r > r esc , where r esc is the solution of the equation 
e esc( r ) = e - Let us now assume for simplicity, that photons which form the 'new' peaks originate 
just at t — r esc (e) (this is a reasonable assumption, especially for the models with no acceleration). 
With the coordinate r of such an emitting ring determined, we now calculate the corresponding 
opening angle (for dipolar magnetic field lines) and then A peak . In Fig.l the blue lines represent 
the values of A pcak (e) found semi-analytically, while the red dots are the Monte Carlo results. 
This branch of solution intersects the horizontal line set by A peak = 0.43 at et urn — 0.9, 4.5, and 3 
GeV for models A, B, and C, respectively. 

For model D, with uniform distribution of primary electrons over the polar cap (but otherwise 
identical to model A) changes of A peak above eturn , occur in the opposite sense. Unlike in previous 
models, here both peaks of the pulse are formed by photons which were emitted along magnetic 
field lines attached to the polar cap at some opening angle 0i n it < 9 pc . These photons are less 
attenuated than those coming from the outer rim, and in consequence the peak separation drops. 
Similar behaviour was obtained by Miyazaki & Takahara (1997) in their model of a homogeneous 
polar-cap. 
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Regardless the actual shape of the active part (i.e. 'covered' with primary electrons) of the 
polar cap (either an outer rim, or an entire cap, or a ring), one does expect in general strong 
changes in the peak separation to occur at photon energies close to high-energy spectral cutoff due 
to magnetic absorption. To illustrate this, the lower panels of Fig. 1 present the energy output per 
logarithmic energy bandwidth at the first peak as a function of photon energy e for models A, B, 
C and D. 



4. Summary 

Motivated by the recent suggestion (Kanbach 1999) that the peak separation in gamma-ray 
pulses in the Vela pulsar may be energy dependent over the range of 50 MeV to 9 GeV, we 
calculated gamma-ray pulses expected in polar-cap models with magnetospheric activity induced 
by curvature radiation of beam particles. Two types of geometry of a magnetospheric column 
above the polar cap were assumed: a hollow-cone column attached to the outer rim of the polar 
cap and a filled column. Four models were considered with three scenarios for beam acceleration. 
The emission showing up as double-peak pulses is a superposition of curvature radiation due to 
beam particles and synchrotron radiation due to secondary pairs in cascades. The changes 
in the peak separation were investigated with Monte Carlo numerical simulations and then 
reproduced (to some extent) with semi-analytical methods. 

We found that for the energy range e < e tU rn ~ a few GeV the peak separation A peak either 
slightly decreases with increasing photon energy e at the rate consistent with Kanbach (1999), 
or stays at a constant level. The gamma-ray emission in this range is dominated by synchrotron 
radiation in all four models considered. The actual behaviour of A pcak depends on physical 
properties of the pairs as well as on their spatial extent in the magnetosphere, which vary from 
one model to the other. 

At e ~ eturn the peak separation A peak makes an abrupt turn, and for e > e tU rn it changes 
dramatically. It increases in the hollow-cone models (A, B, C) and decreases in the filled-column 
model (D), at a rate ~ 0.28 phase per decade of photon energy. This is due to magnetic absorption 
effects (jB — ► e^ 1 ). The numerical behaviour of A pcak in the hollow-cone models was reproduced 
with good accuracy with a simple semi-analytical approach to the condition of magnetospheric 
transparency for a photon of energy e, originating at a given point in the dipolar magnetic field 
and propagating outwards at a given direction. 

The value of e tU rn is model-dependent and for the cases considered here it stays between 
~ 0.9 GeV and ~ 4.5 GeV. To find such hypothetical turnover of A peak in real observational 
data would require, however, high-sensitivity detectors, since for e > e tU rn the expected flux 
of gamma-rays drops significantly (Fig. 1, lower panels). If observed, this turnover would be a 
signature of polar cap activity in gamma-ray pulsars, with high-energy cutoffs in their spectra due 
to magnetic absorption. A detailed account of the effects presented above will be given elsewhere 



- 7- 



(Dyks k, Rudak 1999, in preparation) 
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